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Abstract. Recent measurements of cosmic rays by various experiments have found that the 
energy spectrum of cosmic rays is harder in the TeV region than at GeV energies. The origin 
of the spectral hardening is not clearly understood. In this paper, we discuss the possibility 
that the spectral hardening might be due to the effect of re-acceleration of cosmic rays by weak 
shocks associated with old supernova remnants in the Galaxy. 


1. Introduction 

The recent measurements of cosmic rays by the ATIC [1], CREAM [2], and PA ME T; A [3j 
experiments have revealed an anomaly in the cosmic-ray energy spectrum at GeV-TeV energies. 
The measured spectrum is found to be harder in the TeV region than at GeV energies. The 
spectral anomaly is difficult to explain using standard models of acceleration of cosmic rays and 
the nature of cosmic-ray propagation in the Galaxy, which predict a single power-law spectrum 
over the whole energy range. 

Various explanations for the spectral anomaly have been proposed. These include hardening 
in the cosnric-ray source spectrum at high energies @1 10 El, changes in the cosmic-ray 
propagation properties in the Galaxy 01, and the effect of nearby sources unununum ln 
this paper, we discuss the possibility that the anomaly can be due to the effect of re-acceleration 
of cosmic rays by weak shocks during their propagation through the Galaxy mm- The 
model assumes that cosmic rays are accelerated by strong supernova remnant shock waves. 
After acceleration, cosmic rays escape from the remnants and undergo diffusive propagation in 
the Galaxy. During the propagation, some fraction of the cosmic rays may again encounter 
expanding supernova remnant shock waves and get re-accelerated USE]. This encounter is 
expected to be more often with older remnants as they occupy a larger volume in the Galaxy 
as compared to the younger remnants. Therefore, this process of re-acceleration is expected to 
be produced mainly by weaker shocks, and as weaker shocks accelerate particles with a softer 
spectrum, the re-accelerated particles will have an energy spectrum which is steeper than the 
spectrum originally produced by the strong shocks. The re-accelerated component can dominate 
the GeV energy region, while the non-reaccelerated component (hereafter “normal component”) 
dominates at higher energies. In the following, we will show that this effect can explain the 
observed spectral anomaly. 


Re-acceleration can also be produced by the same magnetic turbulence responsible for the 
scattering and spatial diffusion of cosmic rays in the Galaxy. This process, commonly known 
as distributed re-acceleration, have been predicted to produce strong features on some of the 
observed properties of cosmic rays at low energies, such as the peak in the secondary-to-primary 
ratios at ~ 1 GeV/nucleon. However, earlier studies have shown that this kind of re-acceleration 
is not sufficient enough to produce any noticeable features in the cosmic-ray spectrum at high 
energies [TS]. The efficiency of distributed re-acceleration decreases with energy, and its effect 
becomes negligible at energies above ~ 20 GeV/nucleon. On the other hand, for the re¬ 
acceleration by old supernova remnants shock waves discussed in this paper, the efficiency does 
not depend strongly on energy. It depends primarily on the rate of supernova explosions in the 
Galaxy and the volume occupied by the supernova remnants. Its effect can be larger at higher 
energies as compared to the effect of distributed re-acceleration. 

2. Cosmic-ray transport equation with re-acceleration 

The re-acceleration of cosmic rays by weak shocks in the Galaxy is incorporated in the cosmic- 
ray transport equation as an additional source term following a power-law momentum spectrum 
jjl9j. The cosmic-ray transport equation including diffusion, re-acceleration and interaction losses 
then follows, 
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where a cylindrical spatial coordinate (r, z) system is adopted with the origin located at 
the center of the Galaxy. N(r,p) represents the differential number density of nuclei with 
momentum/nucleon between p and p+dp, D(p) is the diffusion coefficient, and Q(r,p) represents 
the source term. The first term on the left-hand side of Eq. (1) represents diffusion. The 
second and third terms represent losses due to the inelastic interaction of cosmic rayss with 
the interstellar matter and due to the re-acceleration to higher energies respectively, where h 
represents the averaged surface density of interstellar atoms, v(p) the particle velocity, a{jp) 
the inelastic collision cross-section, and £ corresponds to the rate of re-acceleration. The fourth 
term containing an integral represents the population of particles produced by the re-acceleration 
of low-energy particles, assuming that particles are instantaneously re-accelerated to generate 
a power-law distribution of index s. The effects of ionization losses and convection due to 
the Galactic wind are neglected as these effects are mostly important at energies below ~ 1 
GeV/nucleon. Our calculation concentrates only at energies above 1 GeV/nucleon. 

The propagation region is assumed to be bounded at z = ±H, and unbounded in the r 
direction. The matter and the sources are assumed to follow a uniform distribution in the 
Galactic disk with a radial size R. For cosmic-ray primaries, the source term in Eq. (1) is taken 
as Q(r,p) = vO[R — r]H[p — po\Q(p), where v represents the supernova explosion (SNe) rate 
per unit surface area of the disk, H(m) = 1(0) for rn > 0(< 0) is the Heaviside step function, 
and po is the lower cut-off in momentum/nucleon that has been introduced to account for the 
losses due to the ionization process. The source spectrum is assumed to follow a power-law in 
total momentum with an exponential cut-off at high energy. It can be written as function of 
momentum/nucleon as, 

Q(p) = AQ 0 (Ap)- g ex.p ( 2 ) 

where A and Z respectively represents the mass and charge numbers of the nuclei, Q o is a 
constant corresponds to the fraction / of kinetic energy of a supernova explosion that is injected 
into a given cosmic-ray species, q represents the source spectral index, and p c is the cut-off 
momentum for protons at high energies. We have assumed that the maximum total momentum 





for a cosmic-ray nuclei is Z times that of the protons. The cosmic-ray diffusion coefficient 
is assumed to follow D(p ) = Dq/3(p/ po) a , where p = Apc/Ze is the particle rigidity with e 
representing the charge of an electron and c the velocity of light, Dq is the diffusion constant, a 
is the diffusion index, /3 = v/c, and po is a constant. In Eq. (1), the re-acceleration parameter 
is written as £ = rjVv, where V = 47r5ft 3 /3 is the volume occupied by a supernova remnant 
of radius 5ft, and rj is a parameter introduced to take care of the unknown actual volume of 
supernova remnants that re-accelerate cosmic rays. We take 5ft = 100 pc which is approximately 
the size of a supernova remnant of age 10 5 yr expanding in the interstellar medium with an 
initial shock velocity of 10 9 cm s -1 . 

Eq. (1) is solved using the standard Green’s function technique. The solution at r = 0 is 
given by, 
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where Ji is a Bessel function of order 1, 
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and, the function A is given by 
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The local cosmic-ray density is obtained by taking z = 0 in Eq. (3). The first term on 
the right-hand side of Eq. (5) represents the normal cosmic-ray component, and the second 
term represents the re-accelerated component. Re-acceleration takes away particles from the 
low-energy part of the spectrum and puts them into the higher energy region. Therefore, for 
re-acceleration produced by weak shocks where s > q, the re-accelerated component might 
produce visible signatures, such as a bump or enhancement, in the energy spectrum. For re¬ 
acceleration by strong shocks which produces a harder particle spectrum, for instance s = q, 
the re-acceleration effect might not become visible since both the components follow the same 
spectra [19] . 

The spectrum for cosmic-ray secondaries in the Galaxy N 2 (r,p) is obtained by solving a 
similar transport equation to their primaries with the source term given by, 


Q 2 (r,p) = nvi(p)(7i2{p)H[R - r]H\p -p 0 ]IVi(r, p)S(z) (7) 

where vi represents the velocity of the primary nuclei, <j \2 represents the total fragmentation 
cross-section of the primary to the secondary, and N\ is the density of the primary nuclei given 
by Eq. (3). The subscripts 1 and 2 denote the primary and secondary nuclei respectively. 

The secondary-to-primary ratio is calculated by taking the ratio iV 2 / N\. For no re¬ 
acceleration, i.e. £ = 0, it can be checked that Eq. (3) reduces to the standard solution of 
pure-diffusion equation, and the secondary-to-primary ratio becomes inversely proportional to 
the diffusion coefficient at high energies (see e.g., ED]). It can be mentioned that the energy 
spectrum of the secondary nuclei is more sensitive to the effect of re-acceleration than the 
spectrum of the primary nuclei m- Consequently, the secondary-to-primary ratios will also be 
more sensitive to the re-acceleration effect than the primary energy spectrum. For a steeper 
re-acceleration index s > q, the ratio is expected to exhibit an enhancement at lower energies, 
and for a harder index s = q, the ratio will have a flattening at high energies DUES]. For the 
present calculation, since we are considering re-acceleration mainly by old supernova remnants, 
we will only consider the case of s > q with s > 4, which corresponds to a shock Mach number 
of ~ 1.7. 
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Figure 1 . Boron-to-Carbon (B/C) ratio. Solid line: Present work including re-acceleration. 
Dashed line: Pure diffusion model without re-acceleration HU. 
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3. Results and discussions 

We take H = 5 kpc, and the averaged surface density of the atomic hydrogen as n = 7.24 x 10 20 
atoms cm -2 DU- We assume that the interstellar medium consists of 10% helium. The inelastic 
interaction cross-sections are taken to be the same used in the calculations in Ref. El- We take 
the size of the source distribution R = 20 kpc, and the proton low and high-momentum cut-offs 
as po = 100 MeV/c and p c = 1 PeV/c respectively. The supernova explosion rate is taken 
as v = 25 SNe Myr” 1 kpc” 1 , which corresponds to ~ 3 SNe per century in the Galaxy. The 
followings are taken as model parameters: the cosmic-ray propagation parameters (Do,po,a), 
the re-acceleration parameters (?y, s) and the source parameters ( q,f ). 

The values of (Dq, po,a,r],s) are determined based on the measurements of boron-to-carbon 
ratio, and the measured energy spectra of carbon, oxygen, and boron nuclei. The values that 
gives the over-all best agreement between the model calculation and the measurements are found 
to be Do = 9 x 10 28 cm 2 s _1 , p = 3 GV, a = 0.33, ?? = 1.02, s = 4.5. These values correspond to 
the maximum amount of re-acceleration permitted by the available boron-to-carbon data. The 
result on the boron-to-carbon is shown in Figure 1 (solid line), where we have also shown the 
result for the case of pure diffusion (dashed line) with no re-acceleration (rj = 0) taken from Ref. 
HU. The best-fit carbon and oxygen source parameters are found to be qc = 2.24, fc = 0.024%, 
and qo = 2.26, fo = 0.025% respectively, where the /’s are given in units of 10 51 ergs. The 
calculation assumes a force-field solar modulation parameter of cj> = 450 MV. 

The spectra for the protons and helium nuclei calculated using the best-fit values of 
(Dq, po, a,p, s) are shown in Figure 2. The top panel represents proton and the bottom panel 
represents helium nuclei. The data are the same as used in Ref. HU- The source parameters 
required to explain the measurements are found to be q v = 2.21, f p = 6.95% for protons, and 
que = 2.18, fue = 0.79% for helium. It can be seen that the model reproduces the measured 
data quite well, explaining the observed spectral anomaly between the GeV and TeV energy 
regions. 

In Figure 2, the effect of re-acceleration is stronger for protons than for helium nuclei. This 
is partly due to the steeper source spectrum for protons, as compared to the helium nuclei, 
required to reproduce the overall observed spectrum. For a steeper source spectrum, a larger 
fraction of low-energy particles becomes available for re-acceleleration leading into an increase 












Figure 2. Top: Proton spectrum. Bottom: Helium spectrum. The lines represent our results. 
For the data, see the experiments listed in Ref. ttB- 


in the re-accelerated component relative to the normal cosnric-ray component. But, the main 
reason for the stronger re-acceleration effect in the case of protons is due to their smaller inelastic 
collision cross-sections as compared to the helium nuclei. This means that low-energy protons 
can be efficiently re-accelerated to higher energies before they undergo inelastic collision with 
the interstellar matter during their residence time in the Galaxy. For helium nuclei, the loss due 
to the inelastic collision becomes important relative to the re-acceleration. For heavier nuclei 
such as iron for which the inelastic collision cross-sections are even larger, the collision losses 
becomes more important and the re-acceleration effect is expected to be negligible. This is 
shown in Figure 3 where we have compared our model prediction for the iron nuclei with the 
measurements. The calculation assumes qp e = 2.28, and fp e = 4.9 x 10~ 3 %. As expected, the 
re-acceleration effect is hard to notice in Figure 3, and the model prediction above ~ 20 GeV/n 
follows approximately a single power-law unlike the proton and helium spectra. 











Figure 3. Iron spectrum. The line represents our result. For the data, see the experiments 
listed in Ref. mi- 


4. Conclusions 

We have presented a model in which cosmic rays, after acceleration by strong supernova remnant 
shock waves, are re-accelerated by weak shocks associated with old supernova remnants during 
propagation in the Galaxy. The re-acceleration produces a cosmic-ray component that dominates 
the low energy part of the spectrum, while keeping the high-energy part dominated by the 
population that have not suffered re-acceleration. The model is found to explain the cosmic-ray 
spectral anomaly at GeV-TeV energies recently observed for the protons and helium nuclei. The 
study shows that the re-acceleration effect is important only for light nuclei, and is negligible 
for heavy nuclei like iron. Our prediction can be tested by sensitive measurements of heavy 
cosnric-ray species in future. 
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